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ABSTRACT 

Further  studies  of  the  underwater  noise  produced  by  rainfall.  Paul  A.  Elmore,  Hugh  C. 
Pumphrey,  and  Lawrence  A.  Crum;  National  Center  for  Physical  Acoustics;  University  of 
Mississippi;  University,  MS  38677 

A  study  of  the  sound  produced  by  water  drops  striking  a  water  surface  has 
confirmed  some  earlier  results  [Pumphrey  et  al.  J.  Acoust.  Soc.  Am.  85.  1518  (1989)]. 
In  particular,  for  a  certain  well-defined  range  of  drop  sizes  and  impact  velocities,  drops  will 
predictably  and  repeatedly  entrain  bubbles;  this  phenomenon  has  been  named  regular 
entrainment.  In  the  present  study,  various  fixed  drop  diameters  have  been  used  to 
investigate  how  bubble  frequency,  dipole  strength,  and  time  between  drop  impact  and 
bubble  formation  vary  with  impact  velocity.  It  is  found  that  as  impact  velocity  is  increased 
through  the  point  where  entrainment  begins,  both  frequency  and  dipole  strength  decrease  to 
a  minimum  value  and  then  rise  again  as  the  highest  velocity  at  which  entrainment  occurs  is 
approached.  Both  terms  show  increased  variability  near  the  critical  upper  and  lower 
velocities.  The  frequency  tends  to  increase  monotonically  as  drop  size  is  reduced;  drops 
that  entrain  bubbles  at  their  terminal  velocities  tend  to  produce  frequencies  near  14  kHz, 
which  is  also  the  peak  frequency  of  the  natural  rainfall  spectrum.  Finally,  the  time  between 
drop  impact  and  bubble  formation  was  found  to  increase  monotonically  as  drop  size  or 
impact  velocity  increases. 


I 


INTRODUCTION 


This  study  further  investigates  the  underwater  noise  of  rainfall  and  the  phenomenon 
known  as  regular  entrainment.  Previous  studies  (Pumphrev,  Crum,  and  Bjpmp,  1989a; 
Pumphrey  and  Crum,  1989b)  have  shown  that  the  impact  of  a  single  drop  on  a  plane 
surface  of  water  can  make  two  sounds:  first,  an  impact  noise  and  then  a  bubble  noise, 
which  originates  from  the  oscillations  of  a  bubble  entrained  underneath  the  surface  by  a 
drop.  Therefore,  the  spectrum  of  underwater  rain  noise,  which  is  shown  in  Fig.  1,  can 
only  be  caused  by  these  two  effects.  However,  the  impact  noise  is  not  believed  to  be  a 
major  contributor  to  the  total  spectrum.  Although  not  every  drop  entrains  a  bubble,  the 
individual  contribution  by  the  bubbles  are  much  larger  than  that  due  to  the  impact  alone. 
The  majority  of  these  bubbles  are  thought  to  be  oscillating  at  14  kHz  because  the  spectrum 
has  a  large  peak  at  that  frequency. 

These  studies  have  also  shown  that  drops  will  entrain  bubbles  predictably  and 
repeatedly  within  a  well  defined  range  of  impact  velocities  and  drop  sizes,  as  shown  in  Fig. 
2a.  Entrainment  that  occurs  in  this  manner  is  called  "regular  entrainment."  Also  shown  in 
this  figure  are  data  from  Scrimger,  which  give  the  sizes  and  number  of  raindrops  measured 
in  three  different  rainfalls  (Scrimger  et  al.,  1987).  The  comparisons  indicate  that  the  drops 
which  regularly  entrain  drops  at  terminal  velocity  are  found  in  abundance  in  natural  rain. 
Thus,  in  order  to  understand  the  mechanism  of  rain  noise,  we  must  first  leam  about  regular 
entrainment. 

The  regular  entrainment  process  starts  with  the  entraining  drop  striking  the  water 
and  forming  a  conically  shaped  cavity.  This  cavity  will  eventually  stop  growing  and 
collapse  after  a  few  milliseconds.  As  it  does,  the  bottom  part  is  pinched  off  from  the  rest  of 
the  cavity,  forming  a  damped  radially  oscillating  bubble  underneath  the  surface.  This 
bubble  has  a  dipole  radiation  pattern  because  it  is  very  close  to  a  pressure  release  surface. 
Another  process,  which  has  been  named  "irregular  entrainment"  by  Pumphrey  and  Crum, 
can  create  bubbles  as  well  but  by  unpredictable  cavity  deformations  underneath  the  water  or 
entrainment  by  drops  formed  from  jets  erupting  from  the  cavity  (Franz,  1959;  Longuet- 
Higgins,  1989;  Pumphrey  and  Crum,  1989b;  Oguz  and  Prospered,  1989). 

The  present  experiment's  purpose  has  been  to  reproduce  these  earlier  findings  and 
to  examine  more  closely  the  role  of  regular  entrainment  in  the  underwater  noise  produced 
by  rainfall.  Thus,  the  experiment  investigated  the  characteristics  of  bubbles  over  the  entire 
entrainment  region.  The  bubbles'  radii,  frequencies,  and  initial  dipole  strengths  were 
measured  along  with  the  amount  of  time  between  drop  impact  and  bubble  formation.  The 
measurements  were  made  by  keeping  the  drop  size  constant  but  varying  the  impact 


velocity.  The  study  was  then  conducted  over  a  wide  range  of  drop  sizes.  Using  the  results 
of  the  study,  we  can  conclude  whether  or  not  rain  noise  is  caused  by  regular  entrainment. 

METHODS 

The  experiments  were  conducted  in  a  water  tank  measuring  a  cubic  meter  in  volume 
filled  with  tap  water  to  a  depth  between  0.85  and  0.9  meters.  The  transducer  system  used 
for  measuring  acoustic  frequency  and  pressure  was  a  miniature  hydrophone  (Briiel  &  Kjaer 
8103)  and  a  charge  amplifier  (Briiel  &  Kjaer  2635).  The  signal  was  then  read  on  a  digital 
oscilloscope  (LeCroy  9400)  which  is  capable  of  analyzing  the  signal  in  various  ways 
(FFTs,  averages,  integration,  etc.).  The  oscilloscope  was  interfaced  by  a  GPIB  to  a 
computer  (DEC  MINC-73)  for  additional  signal  analysis.  A  diagram  of  the  experimental 
set-up  is  shown  in  Fig.  3. 

Drops  of  tap  water  were  dropped  into  the  water  tank  from  a  syringe  suspended 
above  the  center  of  the  tank.  Various  tubings  or  hypodermic  needles  were  attached  to  it  to 
create  different  sizes  of  drops.  Drops  larger  than  1.71  mm  in  diameter  were  produced  from 
allowing  them  to  drip  from  the  end  of  a  needle,  glass  pipet,  or  piece  of  tubing.  Smaller 
drops  had  to  be  produced  by  feeding  the  water  from  the  syringe  to  a  30-gauge  needle 
attached  to  the  cone  of  a  loudspeaker.  The  speaker-needle  system  was  then  vibrated  with  a 
square  wave  signal  while  allowing  water  to  flow  freely  from  the  syringe  to  the  needle's  tip 
(this  method  gave  the  drops  some  but  negligible  amounts  of  horizontal  velocity).  The 
actual  size  of  the  drops,  assuming  that  they  are  spherical,  was  determined  by  measuring  the 
volumes  of  ten  or  fifteen  captured  samples  with  a  microliter  syringe,  convening  the  each 
volume  to  a  corresponding  diameter,  and  then  averaging  the  diameters.  The  standard 
deviation  of  the  measured  diameters  was  found  to  be  quite  small. 

The  impact  velocity,  vi,  of  the  drops  was  determined  by  measuring  the  height  of  the 
fall,  z,  and  calculating  vi  by  using  the  equation 


where  vt  is  the  terminal  velocity.  The  value  of  vt  for  a  given  drop  size  was  found  by 
computing  it  from  a  polynomial  fit  (Dingle  and  Lee,  1972)  of  earlier  experimental  data 
(Gunn  and  Kinser,  1949).  The  polynomial  fit  and  equation  (1)  were  programmed  into  the 
MINC-73  computer  to  ease  the  computation.  The  error  of  impact  velocity  based  on  the 
error  in  measuring  the  height  of  fall  was  not  found  to  be  significant. 


Using  the  impact  velocities,  the  kinetic  energy,  T,  of  each  drop  at  impact  can  be 
found  from  its  classical  form: 


T  =  Jmv? 
2  \ 


(2) 


where  m  is  the  drop’s  mass.  Since  the  drop's  radius,  not  mass,  was  directly  measured,  it 
was  more  convenient  to  express  T  in  terms  of  the  drop’s  radius,  R,  and  vi.  The  drop's 
mass  is  related  to  it's  radius  by  the  relation 


m  =  pV  = 


(3) 


where  p  is  the  water  density,  V  the  drop's  volume.  Substituting  this  expression  for  mass 
into  equation  (2)  gives  the  desired  form, 


T  = 


(4) 


The  frequency  of  the  bubble's  oscillation  was  measured  by  either  of  two  methods. 
One  method  was  to  measure  the  time  between  the  second  and  seventh  positive  voltage 
peaks  (to  make  sure  transient  behavior  was  excluded)  directly  from  the  oscilloscope  screen, 
as  shown  in  Fig.  4,  and  then  multiply  the  reciprocal  of  the  time  by  five.  The  LeCroy  9400 
measures  the  temporal  difference  and  its  reciprocal  between  any  two  points  along  the  trace. 
The  other  method  was  to  interface  the  oscilloscope  to  the  computer  and  use  a  program  that 
displays  numerical  values  of  the  oscillation's  frequency  and  initial  amplitude.  The  latter 
method  was  more  precise  than  the  first,  but  comparisons  between  the  two  showed  no 
significant  difference.  The  damped  signals  produced  usually  had  quality  factors  of  ten  or 
higher. 

From  the:^  e  readings  the  frequency  was  converted  to  a  bubble  radius,  a,  by  the 

relation 


a  = 


V3kRo/p 

~2nf 


(5) 


where  PQ  is  the  static  pressure  of  the  bubble,  p  the  water  density,  and  f  the  bubble's 
frequency.  Since  the  value  of  these  variables  were  known  for  the  environmental  conditions 
during  the  experiment,  the  equation  reduces  to 


a  = 


3.28 
f  . 


(6) 


Measurements  of  the  initial  dipole  strength  began  with  the  measurements  of  the 
initial  amplitude.  These  measurements  were  also  done  by  either  measuring  initial  voltage 
directly  on  the  oscilloscope  screen  or  by  using  the  program.  When  reading  it  directly  from 
the  oscilloscope,  the  voltage  difference  was  measured  between  the  first  negative  voltage 
peak  and  the  second  positive  peak  (again,  to  exclude  any  transient  behavior),  as  shown  in 
Fig.  5.  The  initial  voltage  was  then  converted  to  initial  dipole  strength,  D,  by  solving  for  D 
in  the  formula 


„  _  D  cos  8  e1^1  ■  rfcie-ft1  - r/c) 

p - -  (7) 

where  p  is  acoustic  pressure,  and  0  is  the  angle  between  the  bubble  and  the  hydrophone 
with  respect  to  the  perpendicular  of  the  plane  of  the  water  surface  going  through  the 
hydrophone,  f  is  the  bubble  frequency,  t  the  time  from  bubble  formation,  c  the  velocity  of 
sound  underwater  (~  1500  m/s),  B  the  damping  factor,  and  r  the  distance  from  the 
hydrophone  to  the  bubble.  This  distance  was  set  at  five  centimeters.  Since  0  =  0,  cos  0  = 
1.  Also,  t  -  r/c  *  0,  so  e  ^1  -r/c)  =  1,  and  eif9  =  1  at  a  peak.  Therefore, 

p  =  ^  =  5  Vpp  (8) 


for  the  equipment  that  was  used.  Since  r  =  0.05  m,  the  above  equation  reduces  to 


D  = 


vpp 

4  . 


(9) 


The  temporal  measurement  between  the  drop  impact  and  bubble  formation  was 
facilitated  by  a  function  on  the  oscilloscope  that  averages  the  traces  from  fifty  drop  impacts 
and  subsequent  bubble  noise  and  displays  the  average  as  a  single  trace.  Using  this 
resultant  trace,  the  time  measurement  was  made  between  the  base  of  the  bubble  noise  and  a 
reference  point  in  the  transient  behavior  of  the  bubble's  formation,  as  shown  in  Fig.  6. 


The  dar  analysis  was  done  using  the  Cricket  Graph  software  package  available  for 
Macintosh  ..  mputers. 


EXPLANATION  OF  DATA  SET  1 

The  first  set  of  data  shows  the  final  velocity,  bubble  radius,  frequency,  and  initial 
dipole  strength  for  each  drop  size.  The  drop  diameter,  measured  in  millimeters,  is  listed 
above  each  table.  The  dates  and  times  shown  at  the  top  of  each  sheet  are  not  of  any 
significance.  The  first  column  is  the  drop's  impact  velocity  in  meters  per  second.  The 
columns  following  it  are  the  measured  values  of  the  bubble's  frequency  in  kilohertz,  initial 
dipole  strength  in  newtons  per  meter,  and  radius  in  millimeters.  The  standard  deviations  of 
each  value  follow  it.  Each  row  represents  the  average  of  thirty  trials,  except  when  noted. 

There  are  four  graphs  that  follow  each  data  page.  The  first  three  plot  the 
relationships  between  impact  velocity  and  bubble  size,  frequency,  and  initial  dipole 
strength.  The  fourth  graph  plots  the  initial  dipole  strength  versus  the  frequency. 

DATA  SET  1 


drop  diameter  •  0.79mm 


Tue.  Aug  8.  1989  2:45  PM 


impact  vel. 

I  (kHz) 

sd  o(  f 

D  (N/m) 

sd  of  0  Bubble 

rad/mm 

sd  of  rad 

1 

3.235 

48.350 

16.422 

0.532 

0.185 

0.200 

0.068 

2 

'  -33 

52.088 

21.562 

0.955 

0.330 

0.152 

0.063 

3 

3.231 

58.626 

25.267 

0.730 

0.290 

0.130 

0.056 

4 

3.240 

40.875 

12.S21 

0.723 

0.290 

0.262 

0.080 

f  (KHz)  Bubble  rad/mm 


drop  diameter  =  0.79mm 


7 


0.40 


•  0.35 


0.30 


0.25 


0.20 


0.15 


0.10 


3.232  3.234  3.236  3.238  3.240  3.242  3.244 

impact  vel.  (m/s) 

drop  diameter  =  0.79mm 


•  Bubble  rad/mm 
- term.  vel.  3.243m/s 


f  (KHz) 

term  vel.  3.243  rrv 


impact  vel.  (m/s) 


D  (N/m) 


Thu.  Jul  20,  1989  7:44  PM 


9 


drop  diameter  »  0.87mm 


impact  vel. 

f  (kHz) 

sd  of  1 

D  (N/m) 

sd  of  D  Bubble 

rad/mm 

sd  of  rad 

1 

3.539 

14.331 

0.612 

0.422 

0.070 

0.229 

0.010 

2 

3.526 

16.852 

3.581 

0.318 

0.065 

0.195 

0.041 

3 

3.521 

15.029 

1.197 

0.385 

0.075 

0.218 

0.017 

4 

3.517 

17.973 

16.704 

0.362 

0.080 

0.182 

0.170 

S 

3.51 1 

16.975 

6.032 

0.350 

0.095 

0.193 

0.069 

6 

3.505 

18.329 

6.341 

0.438 

0.092 

0.179 

0.062 

7 

3.498 

16.899 

1 .213 

0.407 

0.058 

0.194 

0.014 

8 

3.489 

18.698 

1 .690 

0.492 

0.107 

0.175 

0.016 

9 

3.479 

21.1 15 

7.920 

0.652 

0.163 

0.155 

0.058 

10 

3.467 

20.579 

2.743 

0.340 

0.060 

0.159 

0.021 

1  1 

3.454 

20.749 

1.315 

0.463 

0.075 

0.158 

0.010 

12 

3.439 

23.081 

2.084 

0.472 

0.075 

0.142 

0.013 

13 

3.421 

19.439 

1.285 

0.365 

0.045 

0.169 

0.011 

14 

3.401 

18.447 

1.460 

0.405 

0.065 

0.178 

0.014 

1  5 

3.373 

27.705 

7.600 

0.873 

0.223 

0.118 

0.032 

1  6 

3.342 

30.355 

5.319 

1.005 

0.170 

0.108 

0.019 

D  (m/s)  D  (m/s) 


drop  diameter  =  0.87mm 


D  (m/s) 

term.  vel.  3.548  m/s 


drop  diameter  =  0.87mm 


1.2 


i.o  H 


0.8  H 


0.6  H 


0.4  H 


•  • 

•  •  • . 


•  D  (m/s) 


0.2 


drop  diameter  ■  0.93mm 


Thu.  Jul  20,  1989  9:12  AM 


impact  vel.  f  (kHz) 


1 

3.103 

49.786 

2 

3.196 

25.592 

3 

3.276 

22.220 

4 

3.343 

16.426 

5 

3.441 

14.933 

6 

3.390 

14.808 

7 

3.485 

13.216 

8 

3.523 

12.866 

9 

3.558 

15.878 

10 

3.583 

12.190 

11 

3.608 

12.662 

12 

3.629 

12.145 

13 

3.647 

11.535 

14 

3.663 

12.143 

1  5 

3.677 

12.584 

1  6 

3.689 

12.167 

17 

3.699 

14.336 

18 

3.708 

12.731 

19 

3.716 

12.021 

20 

3.723 

12.126 

21 

3.729 

12.024 

22 

3.751 

12.156 

sd  of  f  D  (N/m)  sd 


10.003 

0.917 

7.621 

0.980 

8.228 

C.820 

1.734 

0.573 

2.395 

0.417 

3.837 

0.378 

1.205 

0.307 

2.104 

0.388 

18.225 

0.345 

0.948 

0.323 

1.079 

0.378 

0.832 

0.318 

2.548 

0.240 

1.425 

0.193 

2.352 

0.310 

1.168 

0.302 

8.425 

0.315 

1.433 

0.330 

0.456 

0.242 

0.741 

0.200 

0.744 

0.260 

0.154 

0.185 

0  Bubble 

rad/mm 

sd  of  rad 

0.350 

0.066 

0.013 

0.138 

0.128 

0.038 

0.155 

0.148 

0.055 

0.055 

0.200 

0.021 

0.150 

0.220 

0.035 

0.147 

0.222 

0.057 

0.065 

0.248 

0.023 

0.100 

0.255 

0.042 

0.072 

0.207 

0.237 

0.052 

0.269 

0.021 

0.055 

0.259 

0.022 

0.043 

0.270 

0.019 

0.060 

0.284 

0.063 

0.077 

0.270 

0.032 

0.117 

0.261 

0.049 

0.070 

0.270 

0.026 

0.117 

0.229 

0.134 

0.068 

0.258 

0.029 

0.043 

0.273 

0.010 

0.062 

0.270 

0.017 

0.045 

0.273 

0.017 

0.040 

0.270 

0.003 

f  (kHz)  Bubble  rad/mm 


drop  diameter  =  0.93mm 


•  Bubble  rad/mm 
- term.  vel.  3.768  m/s 


drop  diameter  =  0.93mm 


f  (kHz) 

term.  vel.  3.768  m/s 


impact  vel.  m/s 


(W/N)  a  (ui/N)  0 


drop  diameter  =  0.93mm 


impact  vel.  m/s 


drop  diameter  =  0.93mm 


D  (N/m) 


f 


drop  diamotar  »  1.00mm 


Wed,  vlul  19,  1989  5:18  PM 


impact  ve!. 


f  (kHz) 


Sd  of  f  0  (N/m)  sd  of  D  Bubble  rad /mm  sd  of  radius 


1 

2 

3 

4 

5 

6 

7 

8 
9 

1  0 
1  1 
12 
1  3 
1  4 
1  5 


1  6 
1  7 


1  8 
1  9 
20 
21 
22 


3.1  67 
3.27S 
3.368 
3.449 
3.S18 
3.579 
3  531 
3.677 
3.717 
3.753 
3.784 
3.784 
3.811 
3.835 
3.556 
3.874 
3.891 
3.905 
3.918 
3.929 
3.938 
3.983 


20.500 
17.707 
17.866 
1 4. 1 74 
13.000 
12.022 
11.703 
1 1.148 
11.233 
11.401 
11.291 
11.733 
14.365 
11.287 
14.906 
12.321 
12.403 
13.763 
12.023 
12.215 
16.768 
12.217 


5.368 
5.865 
12.250 
3.267 
2.184 
0.282 
0.366 
0.050 
0.063 
0.138 
0.1  53 
0.400 
2.923 
1 .988 
4.876 
1.370 
2.505 
8.163 
0.138 
1.863 
16.502 
0.449 


0.980 

0.500 

0.405 

0.545 

0.435 

0.378 

0.350 

0.307 

0.265 

0.210 

0.021 

0.193 

0.350 

0.240 

0.378 

0.282 

0.253 

0.275 

0.275 

0.193 

0.182 

0.233 


0.190 

0.160 

0.195 

0.185 

0.128 

0.184 

0.165 

0.231 

0.070 

0.252 

0.018 

0.273 

0.028 

0.280 

0.018 

0.294 

0.022 

0.292 

0.013 

0.288 

0.018 

0.290 

0.020 

0.280 

0.109 

0.228 

0.028 

0.291 

0.190 

0.220 

0.065 

0.266 

0.092 

0.264 

0.105 

0.238 

0.058 

0.273 

0.043 

0.269 

0.048 

0.196 

0.040 

0.268 

0.042 

0.061 

0.126 

0.053 

0.042 

0.006 

0.009 

0.001 

0,002 

0.003 

0.004 

0.010 

0.046 

0.051 

0.072 

0.030 

0.053 

0.141 

0.003 

0.041 

0.193 

0.099 
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drop  diameter  =  1.00mm 


•  Bubble  rad/mm 
term.  vel.  4.01 5m/s 


impact  vei.  (m/s) 


drop  diameter  =  1.00mm 


•  f  (kHz) 

term.  vel.  4.01 5m/s 


impact  ve!.  (m/s) 


(w/N)  a  (ui/n)  a 


drop  diameter  ■  1.07mm 


Thu,  Jui  20,  1989  9:41  AM 


vet. 

f  (kHz) 

sd  of  f 

D  (N/m) 

sd  of  D  Bubble 

rad/mm 

sd  of  rad 

2.941 

85.753 

24.074 

1.218 

0.335 

0.038 

0.011 

3.101 

27.208 

14.211 

1.050 

0.193 

0.1  21 

0.063 

3.238 

15.391 

17.903 

0.398 

0.083 

0.213 

0.248 

3.356 

16.186 

18.079 

0.378 

0.080 

0.203 

0.226 

3.459 

12.376 

1.224 

0.318 

0.075 

0.265 

0.026 

3.549 

12.064 

1.409 

0.300 

0.098 

0.272 

0.032 

3.627 

10.455 

1.742 

0.235 

0.020 

0.314 

0.052 

3.696 

11.616 

0.298 

0.320 

0.065 

0.282 

0.007 

3.757 

1 1.080 

0.202 

0.233 

0.062 

0.296 

0.005 

3.810 

1 1.152 

0.101 

0.275 

0.030 

0.294 

0.003 

3.858 

11.764 

0.888 

0.235 

0.050 

0.279 

0.021 

3.900 

15.138 

9.947 

0.372 

0.172 

0.217 

0.142 

3.937 

1 1.346 

0.805 

0.318 

0.048 

0.289 

0.021 

3.971 

14.851 

10.727 

0.280 

0.1  17 

0.221 

0.160 

4.000 

11.506 

1.139 

0.270 

0.050 

0.285 

0.028 

4.027 

18.925 

18.821 

0.275 

0.050 

0.173 

0.172 

4.050 

11.097 

0.359 

0.292 

0.030 

0.296 

0.010 

4.090 

12.949 

2.745 

0.287 

0.138 

0.253 

0.054 

4.071 

11.599 

0.524 

0.195 

0.030 

0.283 

0.013 

4.106 

21.123 

11.572 

0.138 

0.060 

0.155 

0.085 

4.1  21 

32.643 

18.186 

0.117 

0.060 

0.100 

0.056 

4.136 

24.655 

10.258 

0.102 

0.045 

0.133 

0.055 

4.147 

30.209 

14.412 

0.092 

0.025 

0.109 

0.052 

f  (kHz)  Bubble  rad/mm 


drop  diameter  =  1.07mm 


impact  vel.  (m/s) 

drop  diameter  =  1.07mm 


•  Bubble  rad/mm 
term.  vel.  4.249  m/s 


•  f  (kHz) 

*—  term.  vel.  4.249  m/s 


19 


impact  vel.  (m/s) 


(w/N)  a  (ui/N)  a 


drop  diameter  =  1.07mm 


20 


•  D  (N/m) 

- term.  vel.  4.249  m/s 


< 


« 


drop  diameter  =  1.07mm 


D  (N/m) 


10 


[ill] 


Thu.  Jul  20,  1989  10:07  AM 
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drop  diameter  =  1.22mm 


impact  vel 

f  (khz) 

sd  of  f 

D  (N/m) 

sd  of  D  Bubble 

rad/mm 

sd  of  rad 

1 

2.610 

30.372 

12.849 

1.917 

0.149 

0.108 

0.046 

2 

2.71 1 

1  9.664 

9.813 

1.186 

0.230 

0.167 

0.083 

3 

2.804 

15.230 

5.225 

0.922 

0.193 

0.215 

0.074 

4 

2.891 

15.389 

4.861 

0.925 

0.229 

0.213 

0.067 

5 

2.943 

12.266 

2.049 

0.747 

0.158 

0.267 

0.045 

S 

3.022 

11.996 

2.139 

0.461 

0.178 

0.273 

0.049 

7 

3.1  14 

10.558 

1.343 

0.554 

0.124 

0.31 1 

0.040 

8 

3.188 

10.505 

1.695 

0.502 

0.097 

0.312 

0.050 

9 

3.508 

10.312 

0.370 

0.985 

0.166 

0.318 

0.01 1 

1  0 

3.625 

11.857 

1.483 

0.731 

0.158 

0.277 

0.035 

1  1 

3.729 

12.894 

2.298 

0.988 

0.218 

0.254 

0.045 

12 

3.294 

9.814 

1.794 

0.394 

0.063 

0.334 

0.061 

13 

3.436 

10.630 

0.515 

0.769 

0.1  15 

0.309 

0.015 

1  4 

3.562 

10.923 

0.610 

0.599 

0.076 

0.300 

0.017 

1  5 

3.673 

11.363 

1 .547 

0.702 

0.086 

0.289 

0.039 

1  6 

3.772 

11.094 

0.809 

0.470 

0.1  19 

0.296 

0.022 

1  7 

3.817 

12.056 

2.968 

0.823 

0.189 

0.272 

0.067 

18 

3.860 

18.911 

10.219 

0.785 

0.211 

0.173 

0.094 

1  9 

3.900 

32.727 

21.765 

0.635 

0.247 

0.100 

0.067 

20 

3.935 

30.154 

21 .827 

0.740 

0.290 

0.1  09 

0.079 

*  (k,1z)  Bubble  rad/mm 


drop  diameter  =  1.22mm 


impact  vel  (m/s) 


drop  diameter  =  1.22mm 


Bubble  rad/mm 


f  (khz) 


final  vel  (m/s) 


D  (N/m)  D  (N/m) 


drop  diameter  =  1.22mm 


23 


impact  vei  (m/s) 


drop  diameter  =  1.22mm 


2.0 


1.5 


1.0 


f 


0.5 


0.0 


1 


1  0 


100 


D  (N/m) 


D  (N/m) 


f  (khz) 


Thu,  Jul  20,  1989  10:20  AM 


final  vel. 


1 

2.004 

2 

2.058 

3 

2.099 

4 

2.130 

5 

2.179 

6 

2.217 

7 

2.255 

8 

2.292 

9 

2.328 

1  0 

2.363 

1 1 

2.397 

12 

2.431 

13 

2.465 

14 

2.497 

15 

2.529 

1  6 

2.561 

17 

2.591 

18 

2.622 

19 

2.652 

20 

2.695 

21 

2.738 

22 

2.766 

23 

2.794 

24 

2.821 

25 

2.847 

26 

2.874 

27 

2.900 

28 

2.925 

drop  diameter  *  1.71mm 


(  (kHz)  sd  of  f  D  (N/m) 


46.174 

22.574 

1.897 

16.600 

1.505 

1.383 

12.964 

0.705 

1.160 

11.639 

0.538 

1.055 

10.637 

0.418 

0.927 

10.140 

0.423 

0.883 

9.492 

0.365 

0.765 

9.718 

0.526 

0.775 

9.362 

0.340 

0.740 

9.233 

2.850 

0.715 

9.211 

0.317 

0.600 

9.029 

0.169 

0.635 

8.482 

1.582 

0.623 

9.657 

2.164 

0.620 

9.954 

2.579 

0.662 

8.807 

0.199 

0.552 

11.070 

10.737 

0.585 

10.069 

4.126 

0.723 

10.294 

2.728 

0.745 

10.048 

2.780 

0.690 

11.946 

7.300 

0.875 

11.141 

3.031 

0.787 

13.340 

4.817 

0.938 

11.810 

2.631 

0.925 

20.767 

18.640 

1.160 

27.478 

19.067 

1.143 

27.783 

15.139 

1.228 

37.226 

22.462 

1.180 

sd  of  D  Bubble  rad/mm  sd  of  rad 


0.155  0.071  0.035 

0.147  0.198  0.018 

0.117  0.253  0.014 

0.058  0.282  0.013 

0.058  0.308  0.012 

0.080  0.323  0.013 

0.058  0.346  0.013 

0.072  0.338  0.018 

0.065  0.350  0.013 

0.065  0.355  0.110 

0.052  0.356  0.012 

0.048  0.363  0.007 

0.033  0.387  0.072 

0.102  0.340  0.076 

0.122  0.330  0.085 

0.028  0.372  0.008 

0.050  0.296  0.287 

0.152  0.326  0.133 

0.205  0.319  0.084 

0.150  0.326  0.090 

0.480  0.275  0.168 

0.190  0.294  0.080 

0.172  0.246  0.089 

0.225  0.278  0.062 

0.295  0.158  0.142 

0.453  0.119  0.083 

0.417  0.118  0.064 

0.485  0.088  0.053 
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Bubble  rad/mm 


drop  diameter  =  1.71  mm 


(w/n)  a  (uj/n)  q 


Thu,  Jul  20,  1989  10:40  AM 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 
1  5 
16 

17 

18 

19 

20 


drop  diameter  a  2.17mm 

impact  vel.  f  (kHz)  sd  of  f  D  (N/m)  sd  of  0  Bubble  rad/mm  sd  of  rad 


1.716 

24.796 

6.820 

1.794 

12.134 

0.390 

1.845 

9.786 

0.487 

1.893 

8.827 

0.530 

1.940 

9.573 

0.346 

1.986 

8.580 

0.615 

2.031 

8.496 

0.723 

2.074 

8.389 

0.850 

2.117 

8.325 

0.307 

2.158 

8.346 

0.262 

2.199 

8.128 

0.238 

2.238 

8.115 

0.219 

2.277 

8.745 

0.413 

2.315 

9.510 

0.731 

2.352 

8.905 

0.493 

2.389 

9.558 

0.696 

2.424 

11.381 

3.320 

2.459 

12.179 

8.389 

2.494 

19.439 

17.760 

2.528 

15.525 

5.281 

0.510  0.155  0.132 

0.613  0.022  0.270 

0.603  0.033  0.335 

0.647  0.022  0.372 

0.647  0.028  0.343 

0.657  0.025  0.382 

0.640  0.033  0.386 

0.585  0.055  0.391 

0.600  0.052  0.394 

0.713  0.040  0.393 

0.735  0.035  0.404 

0.708  0.022  0.404 

0.743  0.028  0.375 

0.760  0.058  0.345 

0.762  0.037  0.368 

0.818  0.030  0.343 

0.843  0.065  0.288 

0.828  0.1  10  0.269 

0.825  0.230  0.169 

0.973  0.215  0.211 


0.036 

0.009 

0.017 

0.022 

0.012 

0.027 

0.033 

0.040 

0.015 

0.012 

0.012 

0.011 

0.018 

0.027 

0.020 

0.025 

0.084 

0.186 

0.154 

0.072 


Bubble  rad/mm 


0.5  - 


drop  diameter  =  2.17mm 


K 


Bubble  rad/mm 


impact  vel.  (m/s) 


drop  diameter  =  2.17mm 


f  (kHz) 


impact  vel.  (m/s) 


drop  diameter  =  2.17mm 


impact  vel.  (m/s) 


drop  diameter  =  2.17mm 


f  (kHz)  Bubble  rad/mm 


drop  diameter  =  2.23mm 


i 


impact  vel.  (m/s) 


drop  diameter  =  2.23mm 


2.0  2.2  2.4  2.6 


Bubble  rad/mm 


f  (kHz) 


impact  vel.  (m/s) 


<M  CO 


drop  diameter  *  2.57mm  Thu,  Jul  20,  1969  11. '24  AM 


impact  vel. 

f  (khz) 

sd  ol  f 

D  (N/m) 

sd  of  0  Bubble 

rad/mm 

sd  of  rad 

1 

2.208 

9.927 

2.770 

0.950 

0.188 

0.330 

0.092 

2 

2.167 

11.850 

5.573 

1.035 

0.217 

0.277 

0.130 

3 

2.125 

9.586 

1.993 

0.865 

0.217 

0.342 

0.071 

4 

2.082 

9.233 

1 .400 

0.833 

0.200 

0.355 

0.054 

5 

2.038 

8.155 

0.671 

0.713 

0.068 

0.402 

0.033 

6 

1.993 

7.901 

0.354 

0.710 

0.022 

0.415 

0.019 

7 

1.947 

7.202 

0.640 

0.645 

0.075 

0.455 

0.040 

8 

1.899 

7.625 

0.314 

0.720 

0.037 

0.430 

0.018 

9 

1.850 

7.987 

1 .030 

0.765 

0.085 

0.41 1 

0.053 

1  0 

1.800 

8.460 

1.116 

0.840 

0.1 15 

0.388 

0.051 

1  1 

1.747 

9.050 

3.498 

0.973 

0.253 

0.362 

0.140 

12 

1.693 

8.965 

1.546 

0.930 

0.120 

0.366 

0.063 

1  3 

1.637 

11.207 

6.241 

1.028 

0.297 

0.293 

0.163 

1  4 

1.579 

14.563 

3.944 

1.377 

0.333 

0.225 

0.061 

» 
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drop  diameter  =  2.57mm 


N 

SZ 


impact  vel.  (m/s) 

drop  diameter  =  2.57mm 


•  Bubble  rad/mm 


•  f  (khz) 


impact  vel.  (m/s) 


drop  diameter  =  2.57mm 


D  (N/m) 


impact  vel.  (m/s) 


drop  diameter  =  2.57mm 


•  D  (N/m) 


1 


i  o 
1 1 
12 
13 
1  4 
1  5 


drop 

diamotor  «  2.78mm 

Thu,  Jut  20, 

1989  11:47  AM 

final  vet. 

f  (kHz) 

sd  of  f 

D  (N/m) 

sd  of  0 

Bubble  rad/mm 

sd  of  rad 

2.1 10 

9.506 

1.573 

0.000 

0.000 

0.345 

0.057 

2.170 

27.714 

31.057 

0.000 

0.000 

0.118 

0.133 

2.040 

7.595 

0.313 

0.000 

0.000 

0.432 

0.018 

2.000 

8.280 

2.001 

0.000 

0.000 

0.396 

0.096 

1 .950 

8.423 

2.022 

0.000 

0.000 

0.389 

0.093 

1.900 

7.909 

0.409 

0.665 

0.048 

0.415 

0.021 

2.170 

17.070 

6.272 

0.995 

0.362 

0.192 

0.071 

1.850 

7.246 

0.772 

0.632 

0.092 

0.453 

0.048 

1.800 

7.201 

0.674 

0.695 

0.085 

0.455 

0.043 

1.750 

7.126 

1.736 

0.830 

0.092 

0.460 

0.1 12 

1.700 

8.236 

1.085 

0.940 

0.080 

0.398 

0.052 

1 .640 

8.980 

1.400 

0.828 

0.130 

0.365 

0.057 

1.580 

9.871 

2.767 

0.995 

0.223 

0.332 

0.093 

1.520 

16.567 

5.438 

1.510 

0.320 

0.198 

0.065 

1.460 

17.174 

6.586 

1.417 

0.263 

0.191 

0.073 
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NOTE  TO  THE  USE  OF  d=2.78mm 

The  dipole  strength  values,  D,  that  equal  zero  are  artificial;  they  were  inserted 
solely  to  make  the  graphing  work.  For  impact  velocity  =  2.170  m/s,  N  =  75. 


drop  diameter  =  2.78mm 


Bubble  rad/mm 


final  vel.  (m/s) 


drop  diameter  =  2.78mm 


•  f  (Khz) 


1 1 


— i - - - 1 - 1 - 1 - 1 - - — I 

1.6  1.8  2.0  2.2 


final  vel.  (m/s) 


drop  diameter  =  2.78mm 


final  vel.  (m/s) 


drop  diameter  =  2.78mm 


•  D  (N/m) 


drop  diamotor  ■  2.86mm  Thu,  Jul  20,  1939  12:48  PM 


final  vel. 

f  (kHz) 

sd  of  f  Bubble 

rad/mm 

sd  of  rad 

1 

1.620 

7.233 

0.338 

0.453 

0.021 

2 

1.740 

7.968 

0.939 

0.412 

0.049 

3 

1.470 

8.817 

2.520 

0.372 

0.106 

4 

1.860 

12.950 

5.727 

0.253 

0.112 

5 

1.470 

18.217 

9.365 

0.180 

0.093 

6 

1.390 

8.248 

2.943 

0.398 

0.142 

7 

1.520 

6.611 

0.685 

0.496 

0.051 

8 

1.640 

7.174 

0.352 

0.457 

0.022 

9 

1.800 

7.610 

0.322 

0.431 

0.018 

1  0 

1.850 

8.228 

0.940 

0.399 

0.046 

NOTE  TO  THE  USE  OF  D=2.86mm 

Not  all  of  the  data  points  represent  the  average  of  N  =  30  measured  frequencies. 
For  the  two  impact  velocities  =  1.47  m/s,  N  =  20  and  84.  For  impact  velocity  =  1.85  m/s, 
N  =  20. 


f  (kHz)  Bubble  rad/mm 


drop  diameter  =  2.86mm 
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Bubble  rad/mm 


impact  vel.  (m/s) 


drop  diameter  =  2.86mm 


f  (kHz) 


impact  vel.  (m/s) 


(ocasou^^uiv) 


drop  diameter  *  3.07mm 


Thu,  Jul  20,  1989  1:02  PM 


1 


1  0 
1 1 


impact  va).  f  (kHz)  sd  of  f 


1.390 

11.422 

2.271 

1.451 

8.212 

1.883 

1.521 

7.214 

0.422 

1.582 

6.993 

0.620 

1.640 

6.527 

0.420 

1.697 

7.467 

1.425 

1.751 

7.408 

0.233 

1.803 

8.157 

0.484 

1.854 

9.810 

2.311 

1  904 

13.567 

7.235 

1.952 

16.418 

0.594 

D  (N/m)  sd  of  D  Bubble  rad/mm 

1.060  0.215  0.287 

0.855  0.190  0.399 

0.785  0.045  0.455 

0.710  0.085  0.469 

0.565  0.055  0.503 

0.630  0.025  0.439 

0.625  0.030  0.443 

0.625  0.050  0.402 

0.905  0.225  0.334 

0.745  0.205  0.242 

1.140  0.275  0.200 


sd  of  rad 

0.057 
0.092 
0.027 
0.042 
0.032 
0.084 
0.014 
0.024 
0.079 
0.129 
0.1  Co 
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drop  diameter  =  3.07mm 


Bubble  rad/mm 


impact  vel.  (m/s) 


drop  diameter  =  3.07mm 


f  (kHz) 


impact  vel.  (m/s) 


drop  diameter  =  3.07mm 


D  (N/m) 


impact  vei.  (m/s) 


drop  diameter  s  3.07mm 


•  D  (N/m) 


drop  diameter  -  3.34mm 


Thu.  Jul  20.  1989  1:38  PM 
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impact  vel. 

f  (kHz) 

sd  of  f 

D  (N/m) 

1 

1.246 

21.170 

5.790 

1.185 

2 

1.321 

10.536 

1.551 

0.715 

3 

1.391 

7.186 

0.357 

0.710 

4 

1.458 

7.938 

0.630 

0.825 

5 

1.522 

7.985 

0.31 1 

0.820 

6 

1.583 

12.064 

3.968 

0.975 

7 

1.641 

9.687 

0.427 

0.880 

8 

1.698 

17.586 

6.049 

0.512 

sd  of  D  Bubble  rad/mm  sd  of  rad 

0.230  0.155  0.042 

0.045  0.311  0.046 

0.030  0.456  0.023 

0.065  0.413  0.033 

0.035  0.411  0.016  9 

0.135  0.272  0.089 

0.045  0.339  0.015 

0.145  0.187  0.064 


drop  diameter  =  3.34mm 


•  Bubble  rad/mm 


impact  vel.  (m/s) 


drop  diameter  =  3.34mm 


•  f  (kHz) 


impact  vel.  (m/s) 


drop  diameter  =  3.34mm 


impact  vel.  (m/s) 

drop  diameter  =  3.34mm 


cm  r> 


impact  vel.  f  (kHz) 


drop  diameter  -  3.73mm 
sd  of  f  D  (N/m)  sd  of  D 


Thu.  Juf  20,  1989  1:48  PM 
Bubble  rad/mm  sd  of  rad 
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1 

1.207 

9.050 

4.809 

2 

1.289 

6.957 

1.017 

3 

1.391 

6.377 

0.767 

4 

1.459 

6.661 

0.247 

5 

1.522 

7.124 

0.244 

6 

1.583 

13.785 

11.822 

7 

1.642 

15.768 

8.605 

0.797 

0.305 

0.362 

0.193 

0.725 

0.065 

0.471 

0.069 

0.657 

0.110 

0.514 

0.062 

0.672 

0.025 

0.492 

0.018 

0.735 

0.062 

0.460 

0.016 

0.892 

0.195 

0.238 

0.204 

1.440 

0.517 

0.208 

0.114 

drop  diameter  =  3.73mm 


impact  vel.  (m/s) 

drop  diameter  =  3.73mm 


Bubble  rad/mm 


f  (kHz) 


drop  diameter  =  3.73mm 


ft 


i 


» 


> 


ft 


D  (N/m) 


impact  vel.  (m/s) 


drop  diameter  =  3.73mm 


•  D  (N/m) 


CM  CO  '■t 


drop  diameter  ■  3.99mm 


Thu.  Jul  20,  1989  2:04  PM 


impact  vel. 

f  (kHz) 

sd  of  f 

D  (N/m) 

sd  of  D  Bubble 

rad/mm 

sd  of  rad 

1.247 

6.265 

2.088 

0.525 

0.070 

0.524 

0.084 

1.321 

6.761 

0.452 

0.825 

0.080 

0.485 

0.072 

1.392 

6.842 

0.488 

0.835 

0.110 

0.479 

0.070 

1.285 

8.240 

0.448 

0.906 

0.510 

0.398 

0.048 

1.357 

6.085 

0.137 

0.815 

0.025 

0.539 

0.089 

1.426 

6.749 

0.494 

0.815 

0.055 

0.486 

0.072 

1.459 

9.331 

3.556 

0.448 

0.150 

0.352 

0.038 

drop  diameter  =  3.99mm 


1.3  1.4 

impact  vel.  (m/s) 


•  Bubble  rad/mm 


drop  diameter  =  3.99mm 


•  f  (kHz) 


impact  vel.  (m/s) 


drop  diameter  =  3.99mm 


D  (N/m) 


impact  vel.  (m/s) 


drop  diameter  =  3.99mm 


cm  «o  in  to 


drop  diameter  -  4.48  mm 


Thu.  Jul  20.  1989  2:29  PM 


impact  vel. 

f  (kHz) 

sd  of  f 

D  (N/m) 

sd  of  D  Bubble 

rad/mm 

sd  of  rad 

1  1.035 

12.377 

10.077 

0.963 

0.352 

0.265 

0.216 

2  1.124 

8.488 

2.826 

0.845 

0.362 

0.386 

0.129 

3  0.987 

10.790 

3.400 

0.922 

0.295 

0.304 

0.096 

4  1.167 

8.082 

2.680 

0.479 

0.080 

0.406 

0.135 

Q  5  1.081 

8.132 

1 .576 

0.805 

0.165 

0.403 

0.078 

6  1.208 

21.017 

18.564 

0.825 

0.198 

0.156 

0.138 

* 


r 
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f  (kHz)  Bubble  rad/mm 


drop  diameter  =  4.48  mm 
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impact  vel.  (m/s) 


drop  diameter  =  4.48  mm 


impact  vel.  (m/s) 


drop  diameter  =  4.48  mm 


•  D  (N/m) 


•  D  (N/m) 


drop  diameter  «  4.58mm 


Thu,  Jut  20,  1989  2:44  PM 


impact  vei. 

f  (kHz) 

sd  of  f 

D  (N/m) 

sd  of  0  Bubble 

rad/mm 

sd  of  rad 

1 

1.125 

7.456 

1.788 

0.825 

0.175 

0.440 

0.105 

2 

1.167 

7.650 

1  .864 

0.91  5 

0.155 

0.429 

0.104 

3 

1.146 

6.700 

1.128 

0.825 

0.140 

0.490 

0.082 

4 

1.187 

7.415 

1.759 

0.790 

0.075 

0.442 

0.105 

f  (kHz)  Bubble  rad/mm 


drop  diameter  =  4.58mm 


impact  vel.  (m/s) 


1.18 


1.20 


D  (N/m)  D  (N/m) 


ANALYSIS  OF  DATA  SET  1 


The  data  seems  to  divide  itself  into  three  subsets.  The  first  group  encompasses 
drops  1.22  to  3.73  mm  in  diameter,  the  second  for  drops  larger  than  3.73  mm,  and  the 
third  for  drops  smaller  than  1.22  mm. 

For  the  first  group  the  plot  of  impact  velocity  versus  bubble  radii  or  frequency 
resembles  a  concave  "U"  shape.  As  impact  velocity  increases  from  zero  through  the  point 
where  entrainment  begins,  the  bubble  radius  increases  in  size  to  a  maximum  value  and  then 
becomes  smaller  again  as  the  highest  entrainment  velocity  is  reached.  Regular  entrainment 
discontinues  after  this  point.  The  same  sort  of  effect  occurs  with  frequency,  but  since 
bubble  radius  is  inversely  proportional  to  frequency,  the  concavity  is  reversed.  The  graphs 
of  the  initial  dipole  strength  versus  impact  velocity  are  similar  to  those  for  frequency  versus 
impact  velocity  but  with  a  "V"  shaped  rather  than  a '  U"  shaped  pattern.*  The  variability  of 
each  data  point  is  dependent  on  its  impact  velocity,  being  high  if  the  impact  velocity  is  near 
the  upper  or  lower  velocity  limit  and  low  if  the  impact  velocity  is  in  between  the  two  limits. 
The  final  set  of  graphs  in  this  group-initial  dipole  strength  versus  frequency— cannot  be 
analyzed  alone  because  there  are  not  enough  data  points.  However,  some  features  can  be 
determined  if  they  are  compared  with  Fig.  7,  which  plots  200  drops  of  the  same  size  hitting 
the  water  at  a  constant  speed.  These  graphs  all  show  a  monotonic  rise  in  dipole  strength 
with  increased  frequency.  However,  after  a  certain  point  the  relationship  is  destroyed  and 
prediction  of  initial  dipole  strength  based  on  the  knowledge  of  frequency  is  impossible. 

With  the  exception  of  this  last  graph,  the  graphic  features  of  the  next  two  groups  are 
different  from  first.  The  group  consisting  of  drops  larger  than  3.73  mm  has  more  of  an 
isotropic  behavior.  From  this  point  to  the  point  where  drops  become  too  big  to  cause 
regular  entrainment  the  curvature  in  the  graphs  is  destroyed  and  the  variance  is  constant 
among  the  points  (Thus,  one  can  conclude  that  entrainment  is  unstable  in  this  region).  The 
group  made  of  drops  smaller  than  1.22  mm  in  diameter  only  partially  displays 
characteristics  of  the  first  group.  Since  the  drops'  terminal  velocities  are  lower  than  their 
highest  entraining  velocities,  the  graphs  only  show  part  a  "U"  or  "V".  If  the  drops  could 
fall  faster  than  their  terminal  velocities,  the  patterns  would  probably  be  complete. 

Using  the  highest  and  lowest  entrainment  velocities  observed  in  this  study,  a 
second  graph  of  the  range  of  impact  velocities  and  drop  diameters  is  displayed  in  Fig.  8. 
This  range  matches  the  one  done  in  the  earlier  study,  shown  in  Fig.  2a.  The  curves  are  the 
best  fit  curves  of  the  observed  lowest  and  highest  entrainment  velocities.  From  the 

*  The  marring  of  this  pattern  in  the  graphs  for  drop  diameters  =  1.22  mm  and  2.17  mm  is  probably  due  to 
contaminates  in  the  water. 
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knowledge  of  the  impact  velocities  and  drop  radii,  Fig.  8  can  be  modified  to  show  the 
range  of  regular  entrainment  based  on  kinetic  energy  and  drop  diameter,  which  is  graphed 
in  Fig.  9.  Again,  the  curves  are  best  fit  curves  of  the  highest  and  lowest  entraining  kinetic 
energies. 

From  the  data,  a  determination  of  wnether  or  not  regular  entrainment  causes  rain 
noise  can  be  made.  Fig.  10  shows  an  enlargement  of  Fig.  8  in  the  area  where  the  upper 
and  lower  boundaries  of  the  entrainment  intersect  the  terminal  velocity  curve.  Raindrops 
entraining  bubbles  at  terminal  velocity  are  symbolized  by  diamonds;  the  frequencies  of  the 
bubbles  are  listed  beside  them.  Since  most  of  the  bubbles  entrained  by  drops  falling  at 
terminal  velocity  oscillate  around  12  to  14  kHz,  regular  entrainment  must  be  the  cause  of 
the  14  kHz  peak  found  in  the  spectrum  of  natural  rainfall  and  the  source  of  rain  noise 
underwater. 

At  one  point  during  the  study,  the  peak  was  believed  to  be  caused  by  the  14  kHz 
bubbles  being  the  most  in  number  and  the  loudest.  However,  for  bubbles  entrained  by 
drops  falling  at  terminal  velocity,  the  dipole  strengths  of  14  kHz  bubbles  are  usually  lower 
than  those  oscillating  at  20  to  40  kHz.  Therefore,  the  peak  must  be  due  only  to  the  amount 
of  bubbles  attributed  to  it. 

Figures  11  and  12  display  the  data  in  the  form  of  contour  graphs.  The  first  one, 
showing  similar  regions  of  initial  dipole  strength  for  the  entrainment  region,  indicates  that 
dipole  strength  weakens  as  drop  diameter  is  decreased.  The  other  graph,  displaying  similar 
areas  of  frequency  or  bubble  radius,  shows  that  frequency  steadily  increases  as  drop 
diameter  decreases. 

EXPLANATION  OF  DATA  SET  2 

The  second  set  of  data  measures  the  impact  velocity  and  time  required  for  the 
bubble  to  form  after  drop  impact.  Impact  velocities  are  measured  in  meters  per  second  and 
time  is  in  milliseconds.  As  before,  each  row  represents  the  average  of  thirty  trials.  After 
each  data  sheet  are  two  graphs  which  plot  both  the  relationship  between  the  time  and  the 
impact  velocity.  However,  the  top  graph  is  a  logarithmic  curve  fit  to  the  data  while  the 
bottom  is  a  linear  fit.  Again,  the  drop  diameter  is  listed  at  the  top  of  each  table,  and  the 
dates  and  times  are  unimportant. 


DATA  SET  2 


drop  diameter  •  0.79mm 


Tue,  Aug  15,  1989  2:33  PM 


impact  vel.  time  (ms) 


1 

3.240 

4.641 

2 

3.235 

4.620 

3 

3.233 

4.584 

4 

3.231 

4.569 

time  (ms)  from  drop  impact  to  bubble  formation  time  (ms)  from  drop  impact  to  bubble  formation 


drop  diameter  =  0.79mm 

y  -  -21.744  +  8.145 2x  RA2  =  0.914 


impact  vel.  (m/s) 

drop  diameter  =  0.79mm 


y  =  5.5704e-3  *  xA5.721 8  RA2  =  0.91 3 


time  (ms) 

term.  vel.  3.243  m/s 


i 


< 


drop  diameter  =  0.8 7mm 


impact  vel.  time  (ms) 


Fri,  Jut  21,  1989  11:04  AM 
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1 

3.526 

5.490 

2 

3.517 

5.737 

3 

3.512 

5.708 

4 

3.506 

5.740 

5 

3.499 

5.616 

6 

3.490 

5.558 

7 

3.480 

5.462 

8 

3.469 

5.496 

9 

3.455 

5.506 

1  0 

3.439 

5.464 

1  1 

3.421 

5.204 

12 

3.401 

5.263 

13 

3.373 

5.214 

14 

3.342 

5.152 

i  5 

3.306 

5.1  08 

1  6 

3.263 

5.102 

1  7 

3.539 

5.604 

time  (ms)  from  drop  impact  to  bubbl9  formation  time  (ms)  from  drop  impact  to  bubble  formation 


drop  diameter  =  0.87mm 
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y  =  0.79549  *  xA1 .5541  RA2  =  0.831 


term.  vel.  3.548  m/s 
•  time  (ms) 


drop  diameter  =  0.87mm 

y  =  -  3.0285  +  2.4584X  RA2  -  0.826 


impact  vel.  (m/s) 


drop  diameter  •  0.93mm 


Thu,  Jul  20,  1989  10:01  PM 


imapct  vel.  time  (ms) 


1 

3.204 

5.847 

2 

3.253 

6.155 

3 

3.331 

6.216 

4 

3.390 

6.265 

5 

3.485 

6.763 

6 

3.523 

6.679 

7 

3.558 

6.881 

8 

3.583 

7.081 

9 

3.608 

6.720 

1  0 

3.629 

6.769 

1  1 

3.647 

6.804 

1  2 

3.663 

6.991 

13 

3.677 

7.044 

1  4 

3.689' 

6.886 

time  (ms)  nom  drop  impact  to  bubble  formation  tlme  (ms)  from  drop  impact  to  bubble  formation 


drop  diameter  =  0.93mm 

y  =*  1 .8427  *  xA1 .01 57  RA2  =  0.844 


drop  diameter  =  0.93mm 

y  =  9.0771  e-3  +  1 .8778x  RA2  =  0.832 


impact  velocity  (m/s) 


drop  diameter  -  1.00mm 


Fri,  Jui  21,  1989  11:20  AM 


impact  vel.  time  (ms) 


1 

3.167 

5.847 

2 

3.275 

6.453 

3 

3.368 

6.592 

4 

3.449 

6.714 

5 

3.518 

6.768 

6 

3.579 

6.952 

7 

3.631 

7.030 

8 

3.677 

7.692 

9 

3.717 

7.696 

1  0 

3.753 

7.826 

1  1 

3.784 

7.986 

12 

3.811 

8.113 

1  3 

3.835  . 

7.848 

14 

3.856  ' 

8.188 

1  S 

3.874 

8.126 

1  6 

3.905 

8.075 

17 

3.918 

8.062 

18 

3.929 

7.736 

1  9 

3.938 

7.627 

20 

3.983 

7.522 

drop  diameter  =  1.00mm 

y  =  1 .3061  *  xA1 .3294  RA2  »  0.846 


impact  vel.  (m/s) 


drop  diameter  =  1.00mm 

y  =  -  2.1907 +  2.6048X  RA2  -  0.825 


impact  vel.  (m/s) 


drop  diameter  ■  1.07mm 


Fri.  Jul  21.  1989  8:15  AM 


final  vel.  time  (ms) 


1 

3.101 

6.573 

2 

3.238 

6.881 

3 

3.356 

7.176 

4 

3.459 

7.232 

5 

3.549 

7.460 

6 

—  3.627 

7.633 

7 

3.696 

7.257 

8 

3.757 

7.587 

9 

3.810 

7.594 

1  0 

3.858 

7.903 

1  1 

3.900 

8.725 

1  2 

3.937 

8.994 

1  3 

3.971  . 

8.639 

1  4 

4.000  ' 

8.756 

1  S 

4.050 

8.703 

1  6 

4.071 

9.261 

1  7 

4.090 

10.013 

1  8 

4.106 

9.531 

19 

4.121 

9.607 

20 

4.138 

9.567 

21 

4.147 

9.509 

time  (ms)  from  drop  Impact  to  bubble  formation  time  from  drop  Impact  to  bubble  formation 


drop  diameter  -  1.71mm 


Fri,  Jul  21.  1989  11:25  AM 
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impact  vel 

time  (ms) 

1 

2.004 

11.560 

2 

2.058 

11.840 

3 

2.099 

12.000 

4 

2.130 

12.220 

S 

2.179 

12.420 

6 

-"2.217 

12.600 

7 

2.255 

12.760 

8 

2.292 

12.940 

9 

2.328 

13.140 

1  0 

2.363 

13.240 

1  1 

2.397 

13.400 

12 

2.431 

13.560 

13 

2.465 

13.700 

14 

2.497 

14.180 

1  5 

2.529 

14.220 

1  6 

2.561 

14.360 

1  7 
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time  (ms)  from  drop  impact  to  bubble  formation  ,ime  (ms)  from  drop  impact  to  bubble  formation 


drop  diameter  =  1.71mm 

y  =  6.5009  *  xA0.83239  RA2  =  0.992 


drop  diameter  =  1.71mm 

y  =  2.3442  +  4.6332x  RA2  =  0.990 
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drop  diameter  «  2.17mm 


Fri,  Jul  21.  1989  11:27  AM 


impact  vel. 

time  (ms) 

1.794 

14.720 

1.844 
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1.940 

15.780 
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2.031 

16.230 
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drop  diameter  =  2.17mm 


=  9.6971  *  xA0.72782  RA2  =  0.999 


time  (ms) 


drop  diameter  =  2.17mm 

y  =  4.671 1  +  5.691 6x  RA2  =  0.998 
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impact  vel.  (m/s) 


drop  diameter  «  2.27mm 


Fri,  Jul  21,  1989  11:29  AM 


impact  vel.  time  (ms) 
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1  6 

2.392 
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17 

2.428 

18.800 

18 

2.498 

18.980 

I 

I 


drop  diameter  =  2.27mm 

y  =  9.2074  *  xA0.80057  RA2  -  0.996 


Data  from  "drop  diameter  =  2.27mm" 

y  =  3.3340  +  6.3466x  RA2  =  0.996 


impact  vel.  (m/s) 


o  ro 


drop  diameter  «  2.57mm 


Fri,  Jul  21.  1989  11:31  AM 


.mpact  vei,  time  (ms) 


1  1.519 

1.579 
1.637 

4  1.693 

5  1.747 

6  1.800 

7  1.850 

8  1.899 

9  1.947 

10  1.993 

11  2.038 

12  2.125 

13  2.164 

1  4  2.208 


15.580 

16.110 

16.710 

17.160 

17.560 

17.930 

18.270 

18.864 
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19.584 

19.748 

20.424 

20.680 

20.984 
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t.me  (ms)  from  drop  impact  to  bubble  formation  ,ime  (ms)  from  drop  impact  to  bubble  formation 


drop  diameter  =  2.57mm 

y  =  11.272  *xA0.79038  RA2  =  0.998 


time  (ms) 


impact  vel.  (m/s) 

drop  diameter  =  2.57mm 

y  =  3.901 5  +  7.7893X  RA2  =  0.997 


time  (ms) 


impact  vel.  (m/s) 
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drop  diameter  -  2.06mm 


Fri,  Jul  21.  1009  11:36  AM 


impact  vel.  time  (ms) 
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1  1 
12 


1.520 
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1.501 

10.532 

1.639 

10.930 
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20.240 
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21.230 
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2.042 
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22.300 
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time  (ms)  from  drop  impact  to  bubble  formation  time  (ms)  from  drop  impact  to  bubble  formation 


drop  diamotor  =  3.03mm 


Thu.  Ju!  20.  1989  9:08  PM 


impact  vel  time  (ms) 
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1.457 

19.040 

2 
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19.568 

3 
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20.160 

4 
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20.880 

5 

1 .696 

21.256 

6 

—  1.751 

21.784 

7 
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22.180 

8 

1.879 

22.456 

9 
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1  0 
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1  1 
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time  (ms)  from  drop  impact  to  bubble  formation  time  (ms)  from  drop  impact  to  bubble  formation 


CM  r>  rt  V) 


drop  diameter  =  3.34mm 


Fri.  Jut  21,  1989  11:40  AM 


impact  vel. 


time  (ms) 


1.321 

20.032 

1.391 

20.844 

1.458 

21.384 

1.522 

22.036 

1.583 

22.916 

1.641 

23.476 

1.698 

23.968 
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time  (ms)  from  drop  impact  to  bubble  formation  (jme  (mS)  from  drop  impact  to  bubble  formation 


drop  diameter  =  3.13mm 


y  =  14.697  *xA0.67948  RA2  =  0.993 


1.5  1.6  1.7  1.8  1.9  2.0  2.1 


impact  vel.  (m/s) 

drop  diameter  =  3.13mm 

y  =  7.0770  +  8.2283X  RA2  =  0.991 


impact  vel.  (m/s) 


drop  diameter  a  3.43mm 


Thu.  Jul  20.  1989  9:15  PM 
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impact  vel.  time  (ms) 
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4 
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6 
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time  (ms)  from  drop  impact  to  bubble  formation  time  (mS)  from  drop  impact  to  bubble  formation 


drop  diameter  »  3.73mm 


Thu,  Jul  20,  ',989  9:24  PM 
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1.241 
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1.391 
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1 .459 
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—  1.491 
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7 
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24.864 
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9 
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drop  diameter  =  3.73mm 

20.165  *xA0.55673  RA2  =  0.980 


•  time  (ms) 


- 1 - - - 1 - ' - 1 - - - 1 
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impact  vel.  (m/s) 

drop  diameter  =  3.73mm 
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impact  vel.  (m/s) 


drop  diameter  =  3.99mm 


Fri,  Jul  21,  1989  11:42  AM 
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1.266 
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24.200 
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24.700 

4 

1.426 

25.500 
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1.459 

25.700 

time  (ms)  from  drop  impact  to  bubble  formation  time  (mS)  from  drop  impact  to  bubble  formation 


drop  diameter  =  3.99mm 
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y  =  20.421  *  xA0. 60644  RA2  =  0.969 


impact  vel.  (m/s) 


ANALYSIS  OF  DATA  SET  2 


Figure  13  is  a  comprehensive  graph  of  data  set  2.  It  shows  that  the  time  between 
drop  impact  and  bubble  formation  increases  monotonically  as  drop  size  or  impact  velocity 
increases. 

OTHER  OBSERVATIONS 

Figures  2a  and  8  both  show  the  range  of  impact  velocities  and  drop  sizes  for  regular 
entrainment.  However,  there  are  other  activities  occurring  in  and  around  the  region. 

Just  above  the  upper  velocity  bound  for  drops  larger  than  3.73  mm  in  diameter, 
there  is  a  region  of  irregular  entrainment.  As  the  impact  velocity  steadily  increases  from 
the  lowest  entraining  velocity,  entrainment  eventually  ceases  and  a  lull  of  bubble  formation 
is  observed.  However,  if  the  velocity  is  increased  a  little  more,  entrainment  resumes.  At 
first  glance,  the  entrainment  area  appears  to  have  a  forked  end.  But  what  happens  in  these 
two  areas  is  different.  Below  the  lull,  the  water  cavity  is  in  the  shape  of  a  cone,  and  one 
bubble  pinched  off  the  bottom  every  time.  Above  the  lull,  the  cavity  is  spherically  shaped 
with  bubbles  being  formed  in  an  unpredictable  manner.  Many  times  there  are  several 
bubbles  formed  by  one  drop. 

Another  event  occurs  just  below  the  lower  velocity  boundary  for  drops  about  4.5 
mm  in  diameter.  Here,  surface  bubbles  measuring  about  7  mm  diameter  and  oscillating  at 
1.3  to  1.4  kHz,  are  formed  when  the  drop  hits  the  water.  However,  the  event  is  not  a 
regular  occurrence  and  seems  to  be  dependent  on  contaminants  in  the  water. 

Finally,  when  a  bubble  is  entrained  from  a  drop  falling  at  a  velocity  near  the  upper 
or  lower  velocity  boundary,  the  bubble  will  occasionally  exhibit  second  harmonics,  as 
shown,  for  example,  in  Fig.  14.  The  exact  nature  of  these  higher  harmonics  is  not  yet 
clear,  but  the  most  likely  reason  is  that  the  bubble's  radius  is  behaving  like  a  Duffing 
oscillator.  However,  shape  oscillations  may  be  occurring  some  of  the  time  since  frequency 
peaks  not  corresponding  to  harmonics  have  been  observed  in  some  of  the  bubble's  FFTs. 
Their  occurrence,  however,  is  not  common. 

CONCLUSIONS 

In  summary,  the  characteristics  of  bubbles  inside  the  observed  boundaries  of  regular 
entrainment  have  been  explored.  When  the  drop  diameter  is  kept  constant  and  the  changes 
in  the  bubble's  radius,  frequency,  or  initial  dipole  strength  are  observed  as  the  impact 


velocity  is  continuously  increased  from  the  lowest  to  highest  entraining  velocities,  an 
extrema  is  reached  with  the  variance  in  the  data  points  decreasing  the  closer  they  are  to  it  If 
the  drop  is  large,  (d  >  3.79  mm),  then  any  entrainment  occurs  marginally,  and  the  graphs 
take  on  isotropic  behavior.  If  the  drop  entrains  a  bubole  at  terminal  velocity,  only  pan  of 
these  trends  will  exist. 

As  the  drop  diameter  is  increased  from  the  minimum  to  the  maximum  entrainment 
size,  the  frequency  of  the  bubbles'  oscillations  decrease  from  12  to  14  kHz  to  around  6  or  7 
kHz  .  These  results  show  that  regular  entrainment  is  the  cause  of  the  underw  ater  noise 
produced  by  rain. 

Also,  bubbles  have  been  shown  to  take  longer  to  form  after  impact  as  the  drop  size 
or  impact  velocity  is  increased. 

Finally,  additional  phenomena  have  been  observed  in  and  around  the  entmmment 
region:  irregular  entrainment,  surface  bubbles,  and  non-linear  oscillations.  Although,  the 
details  of  regular  entrainment  has  been  extensively  investigated,  these  other  phenomena 
may  prove  to  be  sources  of  fruitful  investigation  for  future  researchers. 
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FIGURE  CAPTIONS 

Figure  1:  The  sound  spectrum  of  natural  rainfall  made  at  various  windspeeds  between  1.3 
to  1.5  meters  per  second  (Scrimger  et  al.  1987). 

Figure  2a:  A  plot  made  in  a  previous  study  (Pumphrey,  Crum,  and  Bjemp  1989)  of  the 

impact  velocities  and  drop  diameters  necessary  for  regular  entrainment  to  occur. 
Anytime  the  drop  diameter  and  impact  velocity  coordinates  are  located  within 
the  area  bounded  by  the  curves,  regular  entrainment  will  occur 

Figure  2b:  The  graph  plots  the  number  of  different  sized  water  drop  measured  earlier 
(Scrimger  et  al.  1987)  in  natural  rainfall  on  three  different  occasions. 

Figure  3:  A  diagram  of  the  experimental  set-up  used  for  this  study. 

Figure  4:  The  arrows  mark  the  reference  points  in  the  bubble  noise  used  to  measure 
frequency. 

Figure  5:  The  reference  points  in  the  bubble  noise  used  to  measure  initial  voltage.  These 
points  are  marked  by  the  tangendal  intersection  of  the  vertical  lines  with  the 
damped  sine  curve. 

Figure  6:  The  arrows  mark  the  reference  points  used  along  the  drop's  noise  to  measure  the 
time  between  drop  impact  and  bubble  formation.  These  points  are  at  the  base  of 
the  initial  impact  noise  and  a  reference  point  in  the  bubble  formation  sound. 

Figure  7 :  A  plot  of  two  times  the  dipole  strength  versus  frequency  for  200  entrained 

bubbles.  These  bubbles  were  entrained  from  drops  2.57  mm  in  diameter 
striking  the  water  at  2.29  m/s. 

Figure  8:  A  plot  made  during  this  study  of  the  impact  velociti  ‘  drop  diameters 
necessary  for  regular  entrainment  to  occur. 

Figure  9:  A  plot  made  during  this  study  of  the  kinetic  energies  and  drop  diameters 
necessary  for  regular  entrainment  to  occur. 

Figure  10:  An  enlargement  of  figure  5  showing  the  highest  measured  velocities  for 
raindrops.  The  numbers  to  the  sides  show  the  average  of  thirty  measured 
frequencies  in  kHz. 

Figure  11:  Contours  of  frequency  or  bubble  radius  are  shown  within  the  entrainment 
region  up  to  a  drop  diameter  of  3.07  mm. 

Figure  12:  Contours  of  initial  dipole  strength  are  shown  within  the  entrainment  region  uo 
to  a  drop  diameter  of  3.07  mm. 

Figure  13:  Some  of  the  graphs  from  data  6et  2  are  diown  plotted  together.  The  drop 
diameters  associated  with  each  curve  are  listed  beside  them. 


Figure  14:  A  marginally  entrained  bubble  displaying  primary  and  secondary  harmonics. 
The  top  sweep  is  the  acoustic  pressure  of  the  oscillation,  while  the  bottom 
sweep  is  an  enlarged  view  of  its  FFT.  The  numbers  above  the  two  peaks  in 
FFT  correspond  to  the  frequencies  at  those  points. 
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drop  diameter  =  2.57mm;  impact  velocity  =  2.29  m/s;  N=200 
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